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We report a systematic study of transport properties of nanosytems with charge density waves. 
We demonstrate, how the presence of density waves modifies the current-voltage characteristics. 
On the other hand hand, we show that the density waves themselves are strongly affected by the 
applied voltage. This self-consistent problem is solved within the formalism of the nonequilibrium 
Green functions. The conventional charge density waves occur only for specific, periodically dis- 
tributed ranges of the voltage. Apart from the low voltage regime, they are incommensurate and 
the corresponding wave vectors decrease discontinuously when the voltage increases. 
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I. INTRODUCTION 

Transport properties of nanosystems are very differ- 
ent from those of macroscopic conductors. In particular, 
nonlinear or even irregular current-voltage characteris- 
tics seem to be intrinsic properties of these systems. Al- 
though, a great number of theoretical and experimental 
results concern the quantum dots, there is an increasing 
interest in the transport properties of nanowires and sin- 
gle moleculesiS^i The main reason for such a tendency 
is their possible application in the future electronic de- 
vices. Theoretical analysis of the transport phenomena 
is difficult due to the coupling between a nanosystem and 
macroscopic leads. As a consequence, Coulomb correla- 
tions are usually taken into account only approximately. 

The spatial confinement, that originates from the ge- 
ometry of the nanosystems, may be responsible for an in- 
homogeneity of the charge distribution. Additionally, one 
may expect that phenomena typical for low dimensional 
correlated systems, e.g., charge density waves (CDW), 
occur as well. The most of the research on the conduc- 
tance of low dimensional systems with CDW correlations 
focus on the sliding density wavesASil This transport 
mechanism sets on for incommensurate CDW when the 
applied voltage exceeds the depining threshold. For a fi- 
nite incommensurate CDW systems it has been shown 
that the transport properties are controlled predomi- 
nantly by the leads & It has also been shown that the 
conductance of a commensurate CDW system is very dif- 
ferent from that of an incommensurate one. The theoret- 
ical description of incommensurate CDW is similar to a 
theory of the Luttinger liquid. The resulting temperature 
dependence of the conductance is much simpler than that 
of commensurate CDW£ These results lead straightfor- 
wardly to a question abouth the mechanisms that de- 
termine the commensurability of CDW in meso- and 
nanoscales, e.g, whether commensurability of the charge 
distribution depends on the bias voltage and the geome- 
try of the nanosystem. Recent self-consistent investiga- 
tions of the molecular chain show that charge distribution 
strongly depends on the applied voltage i 1 ^ 11 This effect 
has already been observed in molecular devices ii£ In 



analogy to this result, one may expect that the applied 
voltage changes also the charge distribution in the CDW 
system, what affects the current-voltage characteristic. 
The current itself can modify the charge distribution as 
well. As a result, we may obtain a system, where small 
changes of the applied voltage can, through the modifica- 
tion of the charge distribution, drive the system between 
insulating and metallic states. It is possible, that this 
tempting feature of the CDW nanosystems could be ap- 
plied in switching devices. 

In the present paper we use the formalism of nonequi- 
librium Green functions to analyze the charge distri- 
bution in a nanosystem coupled to leads. We demon- 
strate that the applied voltage can induce a transition 
between commensurate and incommensurate CDW. We 
analyze how this transition depends on the geometry of 
the nanosystem. The transport properties determined 
for conventional CDW are compared with the results ob- 
tained for systems with unconventional density waves^ 
In particular, we analyze density waves state with d-wave 
symmetry (DDW) that has intensively been investigated 
as a possible scenario of the pseudogap phase in high- 
temperature superconductors^ Such an order has also 
been proposed as a low temperature phase of some quasi- 
two-dimensional organic conductors^ 



II. MODEL AND METHOD 



The Hamiltonian of the system under consideration 
consists of three parts which describe the nanosystem 
itself, macroscopic electrodes and the coupling between 
the electrodes and the nanosystem: H = H nano + H c \ + 
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The electrodes are modeled by a lattice gas 



of noninteracting electrons with a wide energy band: 

#el = Ek,ff,o ( £ k,a - Ma) c t,o;a Ck :^°" wnere A* denotes 

the chemical potential and a G {L,R} indicates the left 
or right electrode, creates an electron with mo- 

mentum k and spin a in the electrode a. At the mean 
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field level the Hamiltonian of the nanosystem is given by 



, r 



(1) 



whereas the matrices G r {ui) and G a (w) consist of re- 
tarded and advanced Keldysh Green functions, respec- 
tively. For the sake of brevity we do not present the com- 
plete set of equations that determines the Green func- 
tions for the mean field Hamiltonian. Instead we refer to 
Ref. 10 for the details. 



where c| CT creates an electron with spin a at site i of 
the nanosystem, = c\ a Ci a and W%j = (— l)*(cf .Cj CT — 

Cj (7 Ci cr )/2. The potentials Ucdw an d ?7ddw describe the 
strength of interactions that are responsible for the for- 
mation of conventional and unconventional charge den- 
sity waves, respectivelyii^ii The coupling between the 
nanosystem and the electrodes is described by: 
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A few remarks concerning the validity of the mean field 
approach are needed at this stage. It is generally believed 
that the mean field approximation is inappropriate for 
low dimensional systems. However, it has been argued 
that this approximation may be applicable at low temper- 
atures for weak interaction and strong coupling between 
nanosystem and electrodes (see Ref. [ijj and the discus- 
sion therein). Therefore, we restrict our considerations 
to the case J7cdw(ddw) < 4£ and assume a relatively 
strong coupling between the electrodes and the nanosys- 
tem. Additionally, the recent investigations of nanor- 
ings with CDW correlations have shown that the mean 
field results qualitatively agree with the exact onesi& (al- 
though the quantitative differences remain significant). 
These results concern the properties of persistent cur- 
rents in a system, where CDW is pinned by impurities. 
In the present case, the coupling to electrodes should play 
a similar role stabilizing CDW. Consequently, we expect 
that the mean-field analysis provides correct qualitative 
results for the transport currents in a CDW system. 

The thermal averages, which occur in the Hamilto- 
nian (JTJ are calculated self-consistently using the lesser 
Keldysh Green functions: 
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After obtaining convergency, one can calculate the cur- 
rent flowing through the nanosystem: 
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du[f h (oj)-f K (uj)} 



xTr r L HG r (w)r R (w)G a (w) 



(4) 



where f a (co) is the Fermi function of electrons in the elec- 
trode a. The elements of the matrix Tl(r) are given by 



(5) 



III. CONVENTIONAL CHARGE DENSITY 
WAVES 

We start our investigations with a one-dimensional 
(ID) nanowire taking into account only the nearest- 
neighbor hopping. The ends of the nanowire are con- 
nected to the macroscopic leads. The only non vanishing 
elements of T's are assumed to be frequency independent 
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= r, where the sites in the 
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chain are enumerated from 1 to N. The difference be- 
tween the lead's potentials gives the voltage applied to 
the nanosystem eV = fit, — /iR. 

In a ID system d-density wave cannot occur and, 
therefore, in this case we restrict ourselves only to the 
conventional charge density wave (ZJddw = 0). Fig. 
shows current -voltage (I — V) characteristics of 20- and 
40-atom chains coupled to the macroscopic leads. Here, 
we compare the characteristics obtained for C/qdw = 1.5t 
with the results for a non interacting system (Uqdw = 0) 

In the uncorrelated case the I — V characteristic con- 
sists of a series of plateaus smoothly connected by steep 
sections. In the following, we refer to these plateaus as 
original plateaus. In the presence of the CDW interaction 
additional smaller plateaus occur instead of these steep 
sections and the resulting I — V characteristic changes 
from a relatively smooth one to a step function. Within 
the original plateaus the characteristic remains almost 
unchanged. These results clearly show that the CDW 
correlations modify the I — V characteristic only for some 
particular values of V. It occurs as a result of a strong 
suppression of the charge density waves by the applied 
voltage in the regions of the original plateaus. In order 
to confirm this statement we have analyzed the spatial 
distribution of electrons over the nanowire. Fig. [21 shows 
the occupation number as a function of position for a 
20-atom chain at various voltages. We have found that 
within the original plateaus electrons are distributed al- 
most uniformly over the system. Contrary to this, in 
the regions where the additional plateaus occur the elec- 
tron density is strongly inhomogeneous and can be de- 
scribed as (in most cases incommensurate) charge den- 
sity waves, i.e, {flic) can be fitted very accurately by 
n° + A cos (Q • Rj + cf>). It turns out that the CDW wave 
vector Q strongly depends on V. For V = the CDW 
is commensurate with the lattice (|Q| = tt). The same 
holds true for a small voltage. However, when the voltage 
increases, the CDW wavelength increases as well. Fig. 
shows how the density wave vector Q and and the den- 
sity wave amplitude A change with the applied voltage. 
Similarly to the I — V characteristic Q(V) is a step func- 
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FIG. 1: Current-voltage characteristics of one dimensional 
chains containing 20 (Figs, la and lb) and 40 (Fig. lc) atoms 
for the temperature of the leads ksT — O.Oli. Figs, la and lc 
correspond to V = O.lt and Fig. lb to F = 0.3t. Continuous 
and dashed lines show results obtained for Ucdw = and 
Ucdw ~ l-5t, respectively. We have denoted 7o = 2et/h and 
Vo = t/e, where t in the nearest-neighbor hopping integrals- 
Arrows labeled as A, B and C indicate voltages for which the 
charge distributions are shown in Fig. [5] 
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FIG. 2: (Color online) Average occupation (n^) obtained for 
a 20-atom chain with C/cdw = 1.5t The voltage is indicated 
explicitly in the figure, whereas the remaining model param- 
eters are the same as in Fig. Results have been fitted by 
the function (rfcj CT ) = n° + A cos (Q • R; + <f>). For clarity of 
the figure the curves are offset. 



obtained for the 40-atom wire. The main difference is re- 
lated to the number and length of the plateaus. Namely, 
the number of the allowed values of Q increases with the 
increasing length of the nanowire. Moreover, the ratio 
of the lengths of the original and CDW-induces plateaus 
decreases when the length of the nanowire increases. It 
suggests that for a sufficiently long nanowire Q should 
become a continuous decreasing function of V. Compar- 
ison of Figs. andinj} illustrates this tendency. In an 
isolated ID system the CDW wave vector can be changed 
by the modification of the Fermi wave vector (|Q| = 2kp), 
that in turn is a sigle-valued function of the occupation 
number. In the present case Q can be changed indepen- 
dently of the concentration of electrons by means of the 
applied voltage. 

In order to get insight into the physical origin of the 
above results, one can consider an isolated chain. In this 
case, the onset of the charge density waves is determined 
by the CDW susceptibility, defined as a retarder equilib- 
rium Green function: 



X (Q,w) = -((A(Q)|At(Q)», 



(6) 



tion, whereas A(V) is approximately a two value periodic 
function. Comparing Figs. ^ and [3] one can see that each 
of the additional plateaus in the I — V characteristic cor- 
responds to a different value of the CDW wave vector. 
The overall behavior of the I — V and Q(V) character- 
istics obtained for the 20-atom chain is similar to those 



where 

and ct creates an electron with momentum k and spin 
cr. In the static case, the CDW susceptibility is given by 
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FIG. 3: Voltage dependence of the density wave amplitude A 
and the wave vector Q for chains containing 20 (upper panel) 
and 40 (lower panel) atoms. Model parameters are the same 
as in Fig. Arrows labeled as A, B and C denote cases 
presented in Fig. |5] 
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FIG. 4: (Color online) CDW susceptibility as a function of 
the chemical potential /i and the wave vector Q. Lighter re- 
gions correspond to larger values of the susceptibility. These 
results have been obtained for an isolated 20-atom chain with 
the nearest neighbor hopping and periodic boundary condi- 
tions. The solid line shows [i = 2£cos(|Q|/2) (see text for 
explanation) . 



Q)= 2_^/(e k+Q )-/(e k ) 



£k+Q - £k 



(8) 



where = 2t cos(|k|) — /i and /(e) is the Fermi distribu- 
tion function. Fig. 0]shows x(Q,^ = 0) as a function of 
the chemical potential fi for a finite system. The maxima 
of the CDW susceptibility occur for such Q's, that both 
the energies in the denominator in Eq.(JHJ) vanish, i.e., for 
£k = and Q = 2k. In the case of an infinite system 
k changes continuously and the first equation has a so- 
lution for an arbitrary value of the chemical potential. 
Then, the maximum of the CDW susceptibility occurs 
for Q that fulfills the condition /i = 2t cos(|Q|/2). On 
the other hand, for a finite system, k takes on discrete 
values and the maxima of the CDW susceptibility oc- 
cur only for specific values of fj,, i.e., when /i is equal to 
one of the energy levels. This feature is responsible for 
the CDW-induced plateaus in the I — V characteristics. 
In a case of noninteracting electrons the steps arise due 
to resonant tunneling through a multilevel quantum sys- 
tem. The steep sections occur when the successive energy 
levels are taking part in the charge transport. However, 
simultaneously the discussed above criterion for the onset 
of CDW is fulfilled. As a result, the CDW gap opens and 
new plateaus occur in the middle of these steep sections. 

But still there is a question concerning the degree of 
steepness of the sections that connect the plateaus. It 
determines the height of peaks in the differential con- 
ductance. The width of the one particle energy-levels 
is related to T. When the system is weakly coupled to 
the electrodes the energy levels are very narrow and the 
I — V characteristic consists of sharp steps, provided the 
temperature is low enough. Increasing of T smooths out 
these steps. It holds true both in the presence and in the 
absence of CDW, what can be inferred from Figs. QJi and 

In the following we investigate the influence of the 
transverse dimension of the nanosystem. For that pur- 
pose we consider a nanowire of a finite width. In such a 
case the description of the coupling between the nanosys- 
tem and the leads becomes nontrivial. We assume a 
simple model in which leads are described by a two- 
dimensional (2D) lattice gas and the hopping between the 
leads and nanosystem is possible only perpendicularly to 
the edge of the nanosystem (see the inset in Fig. 
Then, the nonvanishing elements of the matrices Tl(uj) 
and Fr^u;) can be calculated directly from the Eq. JSJ: 



f a {uj) 



2^ Em 



'S ia S ja COS (k ■ 



Rjj) S(U! - £ k ,a), 



(9) 

where g denotes the hoping amplitude and Rjj = R^ — 
Rj ■ ^iR.(L) is equal to 1 if the site i is located at the right 
(left) edge of the nanosystem and vanishes otherwise. 

Fig. [5] shows the I — V characteristics for a 4-site wide 
and 10-site long nanosystem calculated for Ucbw = 2t 
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FIG. 5: Current-voltage characteristics of a 10 x 4 nanosys- 
tems coupled to the leads, as shown in the inset. Continuous 
and dashed lines show results obtained for Ucdw = and 
Ucdw = 2t, respectively. We have neglected the energy de- 
pendence of T a (it)) and adjusted the value of the hopping 
energy g in such a way that [r a (<*>)] .. = 0.1tSi a (see Eq. 
for the details). 



FIG. 6: Current -voltage characteristics of a 20 x 6 nanosys- 
tems coupled to the leads, as shown in the inset of Fig. 
Continuous and dashed lines show results obtained for 
C/ddw = and C/ddw = 1.6t, respectively. The remaining 
model parameters are the same described in the caption of 
Fig. 



and Uqbw = 0. Contrary to the ID case an important 
difference between the correlated and uncorrelated case is 
visible only for low voltage. We have found that in this 
case there exist a commensurate CDW with the wave 
vector Q = (n,7r). Increasing of V leads to a disappear- 
ance of the CDW ordering. There exist minor differences 
between both the characteristics for larger values of V. 
However, they appear irregularly and are much smaller 
than in the ID case. Therefore, the voltage-induced in- 
commensurate CDW seems to be an intrinsic feature only 
of the ID systems. One may attribute this behavior to 
general properties of the density waves systems. Namely, 
in the 2D case the nesting of the Fermi surface plays a 
crucial role for stability of the CDW phase. The CDW 
wave vector connects the nested parts of the Fermi sur- 
face. For a square lattice with the nearest neighbor hop- 
ping the Fermi surface is perfectly nested only in the 
half-filled case, what corresponds to the commensurate 
Q = (ir, 7r). Other values of Q do not correspond to per- 
fectly nested Fermi surface. Therefore, incommensurate 
CDW in 2D systems is usually less stable than in ID 
cases. We believe that this property is responsible for 
the visibly different I — V characteristics of ID and 2D 
systems. 



IV. UNCONVENTIONAL CHARGE DENSITY 
WAVES 

In the following we extend our analysis by taking 
into account unconventional density waves, where a con- 
densation of electron-hole pairs with non-zero angu- 



lar momentum occurs^ Such a state, with the an- 
gular momentum 1 = 2, has recently been proposed 
as an explanation for the pseudogap phenomena in 
high-temperature superconductors^ Contrary to con- 
ventional density waves the charge is distributed uni- 
formly over the whole system, but there occur orbital 
currents, i.e., the state breaks the time-reversal symme- 
try. This difference may be visible in the transport prop- 
erties because of the interference between the transport 
and orbital currents. Recent developments in the fabri- 
cation techniques allow one to produce nanowires out of 
high-temperature superconductors^ The pseudogap is 
visible in these systems. Moreover, a discrete switching 
noise in the resistance of the nanowires has been observed 
in the pseudogap regime and explained in terms of the 
formation of the stripe phase. It suggests that interest- 
ing phenomena emerge in high-temperature supercon- 
ductors, when one enters the meso- and nanoscales. One 
may also expect that these new results may contribute 
to understanding of the pseudogap phenomenon. In par- 
ticular, a question arises, whether the results presented 
in Ref. |20| can be explained within the unconventional 
density wave scenario of the pseudogap. 

In order to investigate this problem within the for- 
malism introduced above, we have considered a system 
described by the Hamiltonian JI} with Ucdw = 0. Fig- 
ure |S] shows the resulting current- voltage characteristics 
obtained for Uddw = 1.6t and L^ddw = 0. They are 
very similar to those of conventional 2D CDW for low 
and sufficiently high voltages. In the first case, i.e, for 
low voltage the energy gap does not allow for the current 
flow and the system is insulating. In the opposite case, 
the applied voltage destroys both the conventional and 
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FIG. 7: (Color online) Spatial distribution of the orbital cur- 
rents Aj. The voltage is indicated explicitly in the figure, 
whereas the remaining model parameters are the same as in 
Fig. © 



unconventional density waves. However, for the uncon- 
ventional density waves there exists also an intermediate 
regime, where the transport current is finite but the I— V 
characteristic significantly differers from the results ob- 
tained for the uncorrelated system (Uddw = 0). In this 
regime sharp steps in the I — V characteristic occur, in- 
dicating on the rapid changes of the orbital currents dis- 
tribution. Fig. [7| shows the spatial distributions of the 
orbital currents Aj at various voltages, where Aj is given 
by: 

A, = J (W M +* + - Wi, i+S - W it i-f) • (10) 

For a low voltage Aj is almost independent of the lat- 
tice site i, what indicates on a uniform magnitude of 
the orbital currents. The sudden drop of Aj at the sys- 
tem edges originates from the reduced number of the 
neighboring sites [see Eq. (|10fl ]. Increase of V reduces 
the magnitude of the orbital currents. Additionally, this 
quantity becomes spatially modulated, as can be inferred 
from Fig. The modulation is mostly visible in the 
longitudal direction. It is an almost periodic modula- 
tion, with the period decreasing with increasing V. This 
behavior is opposite to the previously discussed ID con- 
ventional CDW, where the period of the charge mod- 
ulation increases with V. There exist lines where Aj 
changes sign, what corresponds to reverted circulation 
of the orbital currents. The voltage induced transitions 
between various distributions of Aj are accompanied by 
sharp steps in the I — V characteristic. Therefore, one 
could speculate that in larger systems such transitions 
may be responsible for the switching noise in the resis- 
tance of the nanowires in the pseudogap phase^S i.e., 
in the phase, that could be described as unconventional 
density waves stated 



V. DISCUSSION AND CONCLUDING 
REMARKS 

In order to investigate the transport properties of 
nanosystems with charge density waves we have applied 
the formalism of nonequilibrium Keldysh Green func- 
tions. Both the conventional and unconventional states 
have been considered. The most of the already published 
results concern the sliding CDW, where it is a priori as- 
sumed whether the density waves are commensurate or 
not. It has previously been shown that the transport 
properties of commensurate and incommensurate CDW 
systems are different. On the other hand, it is known 
that commensurability of an isolated ID CDW system 
depends on the position of the Fermi level, or equiva- 
lently on the occupation number. In the case of transport 
through the nanosystem, its properties are determined by 
the chemical potentials of the left and right electrodes, 
which arc shifted by the applied voltage. We have shown 
that in the case of ID CDW system the commensurability 
changes with the voltage, whereas the average concentra- 
tion of electrons remains unchanged. These CDW states 
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occur only for specific, periodically distributed ranges of 
the voltage. The number of the allowed values of the 
Q vector increases with the increase of the length of 
the system. Therefore, we expect that for a sufficiently 
long nanowire Q should linearly decrease with the ap- 
plied voltage. We have shown that the applied voltage 
affects also the unconventional density waves. In particu- 
lar, it leads to spatial modulation of the orbital currents, 
especially in the longitudal direction. It is a remnant of 
the voltage-dependent charge modulation, that occurs 
in a ID system with conventional CDW. However, con- 
trary to the CDW case the period of this modulation 
decreases with the increase of V. The difference between 
the I—V characteristics obtained for 2D nanosystcm with 
conventional and unconventional CDW is attributed to 
the interference between the transport and orbital cur- 



rents that occurs in the latter case. 

To summarize, we have shown that the properties of 
nanosystems with conventional and unconventional den- 
sity waves strongly depend on the applied voltage. The 
discussed mechanism should be taken into account also 
in the analysis based on the sliding CDW mechanism, 
since the transport properties depend on the commensu- 
rability. 
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